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Abstract 



We introduce quantum walks with a time-dependent coin, and show 
how they include, as a particular case, the generalized quantum walk re- 
cently studied by Wojcik et al. [Phys. Rev. Lett. 93, 180601(2004)] 
■ which exhibits interesting dynamical localization and quasiperiodic dy- 

' namics. Our proposal allows for a much easier implementation of this 

particular rich dynamics than the original one. Moreover, it allows for 
1 an additional control on the walk, which can be used to compensate for 

' phases appearing due to external interactions. To illustrate its feasibility, 

we discuss an example using an optical cavity. We also derive an approx- 
^ 1 imated solution in the continuous limit (long-wavelength approximation) 

which provides physical insight about the process. 

2 ! 1 Introduction 



Quantum walks (QWs) ^ |2] constitute a promising ingredient in the research 
of quantum algorithms jo] but have also an intrinsic interest, reinforced through 
their connection with quantum cellular automata [2] and with phenomena such 
as Anderson localization or quantum chaos @] 03 E] . 

Both in the discrete or continuous version, QWs provide a mean to explore 
all possible paths on a lattice in a parallel way, which is natural for quantum 
evolution, together with constructive quantum interference along the paths. 
Thus they can allow the development of probabilistic algorithms in a more 
efficient way than their classical counterparts |7j. It is therefore crucial to fully 
explore the possibilities offered by QWs, especially in connection with their 
physical implementation. 

Modified QWs can give rise to new physical phenomena, along with more 
efficient algorithmic applications. Different variations of the standard discrete 
time QW have been proposed, including QWs with two entangled particles jH] 
or entangled coins multi-states QWs ^3E]> as well as QWs with alterna- 
tion of different quantum coins in a certain sequence |12| . More interesting to 
us are the generalized QWs that modify the time evolution by the acquisition 



1 



of position-dependent phases by the walker at every step |S] ■ Those gener- 
alizations show phenomena that differ from the typical linear spreading of the 
wave function in the standard QW, such as quantum resonances and dynamic 
localization. 

Within this spirit, we explore a modification on the standard coined QW, 
which consists on the introduction of a time-dependent coin. As we discuss, this 
modification introduces new possibilities on the walk which are worth investi- 
gating. Here we concentrate on a particular time dependent coin that leads to 
QW equations nearly identical to those corresponding to the generalized QW 
introduced by Wojcik et al. (see also 0] Ej). Our approach presents the 
advantage that the corresponding modifications are made on the coin alone, 
which is a simple one-qubit system, in contrast to the original proposal, which 
requires operations to be performed on a large system (the Hilbert space of the 
walking particle). Moreover, using a time-dependent coin can be used as a con- 
trol mechanism to compensate for a phase arising from some external influence. 
We illustrate this idea with an example. 

This paper is organized as follows. In Section 2 we introduce the basic idea 
of a time-dependent coined QW and relate it to previous works. In Section 3, we 
first review the main aspects of the generalized QW introduced by Wojcik et al. 
[Hj, and then show how an extra transformation (with respect to the standard 
coined QW) on the walking particle, can be encoded into a time-dependent coin, 
and show the equivalence between the obtained generalized QW and that of [3| . 
We also discuss the utility of a time-dependent coin as a control mechanism. In 
Section 4, we show how this generalized QW could be implemented in an optical 
cavity. Then, in Section 5, we derive an approximated continuous limit, a long 
wave-length approximation to this time-dependent QW, which is appropriate 
for describing dynamic localization. Finally, in Section 6 we summarize our 
main results. 



2 Time-dependent coined walks 

The standard QW corresponds to the evolution on a one-dimensional lattice of 
a quantum system (the walker) coupled to a bidimensional system (the coin), 
under repeated application of a pair of discrete operators. Let Hp be the Hilbert 
space of the walker, with {|n) , n € Z} a basis of Hp; and let He be the Hilbert 
space of the coin, with basis {\u) , \d)}. The state of the total system belongs to 
the space H = He ® Hp and, at a given time, can be expressed as 

\il> (*)) = ["« (*) l n » u > + d » (*) l n > d >] ■ (!) 

71 

The evolution of the system is governed by two operators: (i) an arbitrary 
unitary transformation C acting on He, which can be any unitary 2x2 matrix 
and is usually chosen as 
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(with p = 1/2 the balanced Hadamard coin H = (a x + a z ) />/2 is recovered); 
and (ii) the conditional displacement operator S acting on Tip 

S \n, u) = \n+ 1, u) , (3) 
S\n,d) = \n-l,d). (4) 

Altogether, they produce the evolution from instant t — 1 to t given by 

\^(t))=SC\^(t-l)). (5) 

In this paper, we introduce the idea of a modified QW, where the coin 
changes during the evolution, i.e. C(t). In this case, the evolution from instant 
t — 1 to t is defined by 

\iP(t)) = SC(t)\iJj(t-l)) (6) 

A particular case of this would be the proposal in |12j , in which two fixed 
standard coins were alternated in a given sequence, leading to a sub-ballistic 
wave-function spreading for some particular choices of the coin series. In order 
to be more specific, we study the effect of a time-dependent coin of the special 
form 

y/T^pe**® -yfpe 
Notice that Q can be obtained as the sequence of two operations, i.e., 



C (t) - ( n X„»*(t) _ /7J P i*(t) ) ■ ( 7 ) 



C(t) = C (t)C, (8) 

with 

( e-^W \ 
Cb(*) = ( o e»(*) ) • 

and C given by (J2J. 

Again, is quite a general function. In this article we shall restrict 

ourselves to a particular case that, as commented above, leads to a generalized 
QW which is nearly identical to that analyzed in Ref. 0. Other possibilities 
will be considered in a future work. 



3 Using time-dependent coins to implement dy- 
namic localization and quasiperiodic dynam- 
ics 

Recently, Wojcik et al. [S] (see also ^ showed that a generalization of the 
QW (GQW in the following) in which a position-dependent phase $ (n) oc n 
was acquired by the walker with each evolution step, produces quasiperiodic 
dynamics and localization effects. There is a physical reason for introducing 
$ (n): the walker is a physical system that evolves in time, and this evolution 
can introduce such phases, via e.g. external interactions. 

Here we show that such generalization can in fact be recast as a QW with a 
time-dependent coin. We concentrate here on the GQW of [JJ], which is equiva- 
lent to that of 0], but our approach can easily be shown to cover also Romanelli's 
et al. proposal ^6 . In fact, the only difference in the dynamical equations, with 
respect to is that the position dependent phase in |S] goes like $ (n) cx n 2 . 
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3.1 GQW 

Let us briefly present the GQW introduced in [Sj. We first define the discrete 
position operator h such that 

n \n) = n \n) , (10) 
and, related to this one, the phase operator 

E = e 100 ™ (11) 
where 4>q is a constant. Following 0, the evolution of the system is governed 

by 

\i>{t)) = SCE \i>(t-l)). (12) 
The state of the system at a given time can be expressed as 

1^ (*)> = [«" (*) i n ' u ) + d "« (*) i n ' d )] ' ( 13 ) 

n 

from which it is easy to obtain 

Un [t) = e H«-D^ (i _ 1) + ^T~p (f - 1)] , (14) 

d n (t) = e^ n+1 ^° u n+ i (t - 1) - VP 4+i (< - 1)] • (15) 

From the solution of these equations one can evaluate the probability of finding 
the walker at the lattice point n at iteration t by using 

V n {t) = \u n {t)\ 2 + \d n {t)\ 2 . (16) 



We now briefly summarize the main features of the solutions of Eqs. (|14I15|1 . 
Wojcik et al. p] found that for rational values of 4>o/2tt dynamical localization, 
shown by a 'quasiperiodic' behavior of the standard deviation a of the probabil- 
ity distribution, is observed during a transient regime, but for long enough times 
a ballistic diffusion occurs. For irrational values of 4>q/2tt, on the contrary, the 
diffusion becomes suppressed, and the walk shows dynamic localization around 
the starting point for arbitrarily long t. Let us consider the case of a rational 
0o/27r in more detail. 

First we notice that the probability distribution V n (t) is invariant under the 
change 

0o^0o + Trfc, with k G Z. (17) 
Then, if we focus on rational values of 4>q/2tt, we can restrict the study to 

<f> = 2tt~, and < - < i (18) 
P P 2 

where q/p is an irreducible fraction. Moreover, the study can be limited to even 
values of p, since given a case q/p with odd p, there is a value (2q — p) / (2p) with 
even denominator leading to the same probability distribution, as a consequence 
of symmetry (|17|) . 

Keeping this in mind, a numerical analysis of Eqs. (|14I15|I shows that, given 
an even p, the solution of the GQW shows a quasiperiod T — p during the above- 
mentioned transient regime. The duration of this transient, i.e., the number of 
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Figure 1: This figure illustrates the localization features of the GQW during the 
transient regime, by showing the probability that the walker returns close to its 
initial position after m quasiperiods. As can be readily observed, the number 
of these quasiperiods before the transient ends increases with p. 

quasiperiods that it exhibits, turns out to be larger the larger is p. In other 
words, the probability that the walker returns to the initial position after one 
quasiperiod, (Pq (mT) ~ 1) with m e N, increases with p, as we show in Fig. 
1. This figure clearly shows that the loss of localization takes place more slowly 
for larger values of p. 

Apart from this oscillation of quasiperiod p, the standard deviation a of the 
probability distribution shows a faster secondary oscillation that depends on 
q. In fact, one finds q secondary oscillations within each period of the main 
oscillation. These secondary oscillations are more pronounced the smaller is 
q and the larger is T. This is clearly appreciated in Fig. 2, where we show 
the evolution of a for T = p = 110 and four different values of q. Notice how 
the GQW returns (only approximately, remind this is a transient behavior) to 
the initial condition, a (t = 0) = 0, when t = mT, oscillating q times between 
t = mT and t = (m+l)T. 

3.2 An alternative approach 

Here we present our alternative approach. Let us define 




(19) 



It is straightforward to show that the following relationship holds 



SEq — CqEqS, 



(20) 



5 




Figure 2 



Figure 2: Here we show how changing the value of q influences the behavior of 
the walker. We have chosen symmetric initial conditions u n (0) = So_ n /^/2 and 
d n {0) = iSo, n /V2, the Hadamard coin (p = 1/2), p = 110, and four different 
values of q. The standard deviation a(t) against time has exactly q peaks within 
one quasipcriod. 
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with S given by Eas.i|3l4fl and E given by Eq. ljTT|) . 

Repeated use of the above expression, together with the evolution Eq. (|12[) . 
leads to a modified form of the evolution equation, which can be expressed as 

i$(t)>= (E y^(t)), (21) 

where \ip (t)) verifies 

|V>(t)> = SC(t)\i>(t-l)), (22) 

|V (0)) = |^(0)>, (23) 

and C'(t) is a time-dependent coin operator, defined as 

^)-«^=(^t ^:r), (24) 



i.e., fcbl = C*o(*) 5 c.f. Eq. J7J), with $(t) = </> i. By this simple procedure we 

have demonstrated that the generalization of the QW introduced in Ref. [H] can 
be obtained by introducing a suitable time dependent coin. Although the proba- 
bility amplitudes are not identical to those of the GQW due to the phase factors 
appearing in Eq. 1)21(1 , the probability distributions are the same obtained either 
with |f/> (i)) or with \ip (t)), and both descriptions are thus equivalent from this 
point of view. 

The equivalence nevertheless breaks down for the QW on the circle, as the 
phase added in [5] depends on the position, so that a difference may arise in the 
circle when passing from position — L to position +L. 

Now we write down explicitly the equations of evolution for our alternative 
approach. By performing a decomposition analogous to Eq. I)13|l . the equations 
of evolution become 



u n {t) = e-^ Vp" n -i (*-!) + ~pd n -x{t-l) , (25) 



d n (i) = e*** [Jl~ P u n+1 (t - 1) - y/pdn+i (t - 1)J . (26) 

which are equivalent to Eas. l|14ll5|) . as Eq. I|21|l provides the connection between 
both descriptions, which reads 

u n (t) = u n (t) e mt ^ , d n (t) = d n (t) e int *° . (27) 

We now transform the coupled equations (|25|) and into space-time re- 
cursive equations for u n and d n , where both components are decoupled. We 
start from 

= S6(t+l)\il>(t)), (28) 

= Ct(t)5t|^(t)>, (29) 

and making use of 

C(t + l)=C C{t), (30) 

we obtain, after some algebra 

ci w (t + 1)> - iv> (t - 1)) = vp E E K-i (*) e ~ lt *° - (*) e ' it0 °] i«. «> 

a—u,d n 

(31) 
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or, equivalently, 

u n (t + l)e^ -u n (t-l) = Vp[ lln _ 1 (t) e - 1 *- Ml (i)e t *],(32) 
d n (t + 1) e"^° - d n (t-l) = y/p [d n -i (t) e~ u *° - d n+1 (t) e«*>] , (33) 

Finally, the probability of finding the walker at the lattice point n at iteration 
t is given by 

V n (t) = \u n (t)\ 2 + \d n (t)\ 2 ee P£(t)+V*(t). (34) 

Since \u n (t)\ = \u n (t)\ and \d n (t)\ = \d n (t)|, there is no difference between the 
probability distribution for the QW on a line calculated with Eqs. 1(251 126(1 or 
with Eqs. ifTHlToj l. as already commented. 

3.3 The time-dependent coin as a control mechanism 

In this subsection we discuss how a time-dependent coin can be used to gain 
control over a possible phase arising during the walk, as a consequence of ad- 
ditional interactions @]. We illustrate this idea with an example which shows 
that, at least in some cases, the position-depending phase acquired between two 
steps in the walk could be canceled by an appropriate action on the coin. 

Let us assume that the walker is subjected to the effect of the GQW defined 
by Eq. 112|) . We have shown in the previous subsection that this kind of 
QW is equivalent (modulo a final phase) to one with a time-dependent coin. 
Intuitively, if one wants to compensate for the phases acquired during the GQW, 
one should replace the coin operator C by a time-dependent operator C(t) 

defined by C(t) = (c'o) In this way, the evolution is governed by 

|V (t))= SC(t)E |V (t-l)). 
Using the properties given in section 3.2, one obtains: 

\1>(t)) = (^'(^y^O)), (35) 

Showing that the combined action of the phase operator Eq and the time- 
dependent coin defined above, is equivalent (up to a phase give by the action of 

(^Et^j ) to the standard quantum walk introduced in section 2. In other words, 
Eq. (|35|l can be written, when decomposed in the \n, u >, |n, d > basis, as 

a n (t) = e int *°a s n (t) 

where a = u,d. The coefficients a n (i) then correspond to Eq. 1)35(1 . whereas 
a s n (t) stand for the standard QW. In this way, the complex dynamics arising 
from the GQW translates into a trivial phase. 

4 Implementing the generalized quantum walk 

Along recent years there have been many proposals for the experimental imple- 
mentation of QWs. These cover both systems whose dynamics can be described 
only within the framework of quantum mechanics El 1161 1171 1181 ITU] as 
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well as setups whose description does not require quantum mechanics |2UI 1211 
[221 EU E] • In fact, the QW on the line was nearly implemented in an optical 
cavity as it was highlighted in [21] and fully discussed in [221 - Although an 
experimental realization of the QW using only classical means has been commu- 
nicated recently , it is a fact that there has been little experimental research 
about this process. 

Here we comment on how the GQW we are studying could be implemented 
in an optical cavity. We follow our approach to the GQW as it is more easily 
implementable than the original proposal by Wojcik et al. [Hj. This is due 
to the fact that with our approach it is only needed to modify the unitary 
transformation acting on the qubit, which is a 2-dimensional system, while the 
original proposal [5] implies acting on all the points of the lattice. 

In |2 11122) . it was shown that the QW on the line can be implemented by the 
frequency of a quasi-monochromatic field, e.g. an optical pulse of appropriate 
duration, inside an optical cavity. As stated, in this classical implementation the 
role of the walker is played by the field frequency, and the role of the coin can be 
played, e.g., by the field polarization. The simplest scheme is that represented 
in Fig. 3 |22i without EOMbis (see below for the role of this element): The elec- 
trooptic modulator (EOM) implements the displacement operator, Eqs. 
by increasing (decreasing) the frequency of the horizontal (vertical) polarization 
component of the field. As for the unitary transformation, C, it is performed 
by a half-wave plate (HWP) with suitably oriented fast axis [27j. Thus in a 
cavity round-trip a step of the QW is performed. The optical cavity allows the 
repetition of the process through feedback. The number of steps of the QW 
that can be implemented depends on factors such as the technical limitations of 
the EOM and on the losses of the cavity (this last factor could be compensated 
by introducing gain in the cavity, as in [23)- We address the interested reader 
to [22] for more details. Let us remark that this simple scheme is very close to 
what was actually performed in the experiment of Bouwmeester et al. |25| (see 
[22] for a full discussion). 

We can take this scheme as a basis for the implementation of the generalized 
QW. In order to perform the GQW, one needs to implement the time dependent 
unitary transformation C (t), Eq. I)24|). This can be done by adding one optical 

element between EOM and HWP to implement 

(C ) , Eq. H9. This 

is the 

role played by EOMbis in Fig. 3. Consider first a single step of the GQW, i.e., 
that corresponding to iteration t. For this t one must implement ypo^j j which 
can be done in a straightforward way: The implementation simply consists 
in the addition (subtraction) of <j>ot to the phase of the horizontal (vertical) 
polarization component of the field. This can be easily carried out, e.g., by 
introducing a second EOM, EOMbis in Fig. 3, to which a suitable (constant) 
voltage is applied. Now, in order to implement C (t), this added (subtracted) 
phase must be increased at each cavity round-trip, what is done by applying a 
staircase voltage to EOMbis (represented in Fig. 3): The voltage must remain 
constant while the light pulse is traversing EOMbis, in order to modify the 
phase and not the field frequency, and then be rapidly increased for the phase 
increment takes the value </>o (t + 1) in the subsequent round-trip. We think 
that this simple scheme, which can be implemented with current technology (it 
consists in adding a single element to the device already used in [213), could 
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Figure 3: The proposed experimental setup for implementing the GQW. See the 
text for the description of the different components. The upper part represents 
the staircase voltage that has to be applied to the EOMbis to implement the 
time-dependent coin, and t is the cavity roundtrip time. 



allow the experimental investigation of the GQW. 



5 A long— wavelength approximation to the gen- 
eralized quantum walk 

Up to now we have shown how the GQW can be alternatively produced by 
means of a time dependent coin, how it could be experimentally implemented, 
and also how the time dependent coin can be used in the presence of phases 
in the walker displacement for controlling or tuning the GQW. Now we will 
try to get some insight into the physics of the GQW, by deriving a continuous 
version through a long-wavelength approximation. In this way we derive a wave 
equation that constitutes a continuous propagation analog of the GQW. The 
analogy helps to visualize the kind of physical process that produces the GQW. 

We introduce here a long wave approximation by following the same lines as 
in j^|2H]. Our starting point is the recurrence equation 



a n (t + 1) - a n (t - 1) = yfp [a n _i (t) e 



-itj> a t 



I'n+l 



(t)e 



(36) 



where a m (t) stands for both u m (t) and d m (t) . In this way, Eq. I|36|l corresponds 
to Eqs. (|32(1 and (|33|l after the factors e ±l ^° on the left hand side have been 
neglected 1 . In |21II28| it was shown that it is necessary to introduce two discrete 

lr This approximation is perfectly justified. Perhaps it is more clearly seen if instead of the 
unitary transformation 1241 . one uses 



C(t) 



V 1 - P e 



-i0 o (t+i) 



-i<Ao(*+|) 



In this case, the exponential factors we are neglecting do not appear on the left-hand side, as 
in Eqs. 1321 1331 . but on the right— hand side of these equations in the form e^ 1 ^ ^^ 2 ) which 
can be approximated by e ±ll ^ ot for large enough t. 
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fields Af; (t) in order to preserve the symmetry of the QW. Thus we define the 
new fields (t) through 



a n (t) = A+ (t) + (-if A' (t) . 
By inserting this definition into Eq. (|36|) . one immediately obtains 

At (t + 1) - At (t - 1) = ±VP [A*-! (t) e-*** - A± + i (*) e^] 
which is convenient to rewrite in the form 



(37) 
(38) 



A±(t + 1)-A±(t-1) = ± y /p[Ati(t)-At 1 (t)}co8ct>Qt 

Ti^P [ A tx (*) + (*)] sin ^- 



(39) 



Denoting by x and t the continuous space and time variables, and by Ax and At 
the spacing between lattice points and time between iterations, we can define 
the adimensional continuous variables £ = x/ Ax and t — tj At and think of 
Eq. (|39|) as the discretization of the following partial differential equation 



E 

k=0 



1 



Q2k+l 



[2k + l)\ dT 2k + X 



1 



d 



2k+l 



= TVPCos(^r)£ (2/t + i)! ^ 2fc+ 



fc=0 



T^si 11 (for) 
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2 A- 



^ (2*)! 



(40) 



which constitutes a continuous limit of the GQW. 

Taking into account Eq. (|37(1 . and the fact that the discrete fields a n (t) 
describe both u n (t) and d n (t), the continuous versions of these fields, which we 
denote by u (£, r) and d (£ , t), are calculated through 



«(£,r) = [/+(e,r) + (-l) i C/-(e,r) ! 
d(£,r) = D+(e,r) + (-!)*£»- (€,t), 



(41) 
(42) 



with [7± (f , t) and L> ± (£, r) the solutions of Eq. gDJ} for A± (£, r) = C/ ± (f , r) 
and A 1 * 1 (£,t) = -D"* 1 (£,t), respectively (see the Appendix). 

The long-wavelength approximation consists in retaining the lowest order in 
Eq. I|40(l . Importantly, we further neglect the third temporal derivative. We 
address the reader to the Appendix for full details. After all of this, we are left 
with 



JV(e,r) = wp 



d 



o 2 



cos (4> t) — + - sin (4> t) — 



on 



1 d 3 



where the new fields B ± (£, t) defined by 
B ± (e,r)= J 4 ± (e,r)exp 



cos(0 o t) 

00 



(43) 



(44) 



have been introduced. This equation can be solved analytically, and the explicit 
solution is derived in the Appendix. 

We notice that Eq. 143|) has time-periodic coefficients and, consequently, as 
we have retained only the first derivative with respect to time, their solutions 
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Figure 4 

Figure 4: A comparison of: (a) an exact numerical evaluation of Eqs. 1(32133(1 
with (b), the long- wavelength approximation derived in this section. Both cal- 
culations correspond to a value <j>o = 27r/150. Only even lattice points have 
been considered. 

are time-periodic. We can then expect that the solutions of Eq. 1(43(1 describe 
approximately the periodic solutions of the GQW, which appear when 0o is 
an irrational multiple of 2n, but not the quasiperiodic solutions (</>o a rational 
multiple of 2tt) except in the cases with very long quasiperiod. Obviously, this 
partial description of the solutions is the price to be paid after neglecting the 
third time derivative in Eq. 1)40(1 . 

Before discussing the physical meaning of Eq. ((43(1 . let us first compare 
the exact solution of the time-dependent coined QW, Eqs. I(32I33|) . with the 
approximated continuous solution we have just derived. In order to do that, we 
have chosen a value for the phase <j>o (tfio — 27r/150) for which the quasiperiod 
T is very large (T = 150 in this case). We have taken symmetrical initial 
conditions too (i.e. uq(0) = 1/V% and do(0) = For the continuous 

version, we will take A ± (^, 0) to be a superposition of gaussians with a width 
w (see the Appendix for details). 

Fig. 4 shows both the exact probability distribution V n (t), Fig. 4(a), and 
the approximated continuous solution, Fig. 4(b), for time running from t = 10 
to t = 140, at intervals of 10 time units. For the exact probability distribution 
only even points of the lattice, for which the probability is non zero, are shown 
and joined for an easier visualization. We have chosen w = 0.65 to evaluate 
V(£,,t). One sees how similar these distributions are, except for t close to T/2, 
where the continuous distribution is wider. Then, the approximated continuous 
solution can be considered as a good approximation for cases with periodic 
behavior or with quasiperiodic behavior with very long transients. 

Fig. 4 is complemented with Fig. 5, where we show: the exact (n,V n (t)) 
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Figure 5: Plots of (n,V n (t)) and (£,"P(£, t)) projections for three different time 
values: t = 20, t = 70, t = 110. 



(with even points joined again) on the top row; the approximated (CP^t)) 
on the bottom row; and finally, in the middle row, the same as in the bottom 
row (i.e., the long-wavelength approximation) but evaluated only at discrete 
position values for a better comparison of the previous two results. We do this 
for three different time values (t = 20, t = 70, t = 110). Again, as in Fig. 4, 
one sees how V n (t) and r) are very similar, except near the semiperiod. 

Finally, we compare the evolution of the quadratic deviation a 2 in position 
using both the exact distribution V n (t) and the continuous distribution V(£,t) 
in Fig. 6. We continue in the dynamic localization case with <f>o = 27r/150, and 
compare the exact case (a) with five continuous limit cases (b) corresponding 
to w = 0.45, w = 0.55, w = 0.65, w = 0.75 and w — 0.85. Notice that the 
behaviors of both the exact case and the continuous limit are similar, except 
for the fact that with the continuous limit one obtains an "excess of quadratic 
deviation" , specially within the proximity of the semiperiod, because of the 
already mentioned problem with the width. 

The above results show that the continuous long-wavelength approximation, 
Eq. (|43|l . is a good qualitative approximation, even a reasonably good quanti- 
tative approximation, to the GQW in the dynamical localization regime. We 
have already commented that the failure in describing the diffusive dynamics 
occurring in the GQW for rational values of 4>o is due to the neglect of the third 
order time derivative in deriving Eq. (|43|l . an approximation made in order to 
obtain analytical expressions. 

Eq. (|43l) is a linear wave equation describing the propagation of waves in 
a medium with special dispersion properties: The dispersion coefficients (those 
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Figure 6 

Figure 6: A plot of the quadratic deviation a 2 in position using both the exact 
distribution V n {t) (left) and the continuous distribution V(£,t) (right), again 
with 0o = 27r/150. On the right panel, the different curves correspond to 
w = 0.45, w = 0.55, w = 0.65, w = 0.75 and w = 0.85. 

multiplying the higher order spatial derivatives) are time periodic, as well as the 
wave group velocity (the coefficient multiplying the first order spatial deriva- 
tive). Thus, the essential for dynamical localization correspond to the vanishing 
of the time-averaged group velocity, while its periodic time dependence is the 
responsible for the "bouncing" of the probability distribution, see Fig. 4. As 
for the rest of spatial derivatives, they introduce a distortion on the probability 
distribution (due to dispersion) that turns out to be reversible again because of 
the time periodicity of the coefficients. Certainly, a group velocity that changes 
its sign periodically is not a common situation for waves, but the analogy that 
Eq. H43fl establishes provides an alternative physical picture that, as we have 
seen, helps to understand dynamical localization in the GQW and could help 
for the search of propagation phenomena in which this phenomenon could mani- 
fest. In this sense, it is interesting to notice the similarity between this equation 
and that describing beam propagation in waveguides with a bent axis |.'i(J|, an 
optical process in which Bloch oscillations and dynamical localization have been 
recently experimentally observed |3J (we note that in |S] the connection between 
the GQW and Bloch oscillations was put forward). 

6 Conclusions 

In this article we have introduced QWs with time dependent coins. We have 
considered a particularly simple case that turns out to be equivalent to the gen- 
eralized QW (GQW) introduced by Wojcik et al. 0. This GQW exhibits very 
striking dynamical properties, particularly dynamical localization. We have 
shown how our alternative proposal can be used as a control mechanism. In 
addition, this time-dependent QW is particularly interesting from the imple- 
mentation point of view, as only simple actions on the coin-qubit are required 
for that. 

We have also obtained a long-wavelength continuous approximation of the 
GQW equations that have allowed us the derivation of an approximated ex- 
plicit continuous solution that works quite well during the dynamic localization 
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regime. The continuous equation from which this solution was derived is a linear 
partial differential equation describing pulse propagation in a dispersive medium 
with periodical time dependence in the dispersion coefficients. This continuous 
limit has lead us to interpret the main feature of GQW, the dynamic localiza- 
tion, as a propagating solution in the dispersive medium with null mean value 
of its group velocity. 
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7 Appendix 



Now we perform the long-wavelength approximation that consists in retaining 
terms up to k = 1 in Eq. (|40|l . i.e. 



d 1 d 3 
Ih + 3! tV3 



Tv^cos^qt) 



Tiy/P sin ((for) 



9 1 9 3 
d£ + 3!9^ 

1+ 2!d£2 



A ± (^t). (45) 



The third-order derivative on the left hand side makes it hard to obtain an 
analytical solution. In the case of the standard QW, the third time deriva- 
tive was approximated by making use of the lowest order expansion (k = in 
Ea. (l4()p ) 2j, but in our case the time-dependent coefficient of the remaining 
linear term renders this approach useless. We then make a further approxima- 
tion and neglect the third order derivative in time. By making the change 



B ± (^r) = A ± (^r)e X p 



00 



one obtains 



9 i d 2 1 d 3 

cos (0 o r) — + - sin (for) ^ + - cos (0 o t) — 



(46) 



B ± &r) 



which we pass to solve here. 

By Fourier transforming Eq. I|45|) one easily gets 



/-t-oo 
dk B ± (fc,0)e lfc «e 
-CO 



(47) 



(48) 



sin {4> t) k cos(0 o r) - l fc2 
0o 200 



^-^U 3 , (49) 



where 



1 /" +00 

B ± (k,0) = — d^B ± (^0) 

2n J-oo 



(50) 



In order to solve the integral 1481) . one must fix B ± (^,0), i.e., A^ (£,0). 
Following [25, we assume that 



A+ (1) - A- (1) ^ A+ (0) + A-(0), 



(51) 
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and then, by using Eq. ((371) . 



^»(0) = -K(0)±a„(l)]. 



(52) 



Notice that a n (l) is evaluated from Eqs. 1)251 126(1 once the initial condition 
a n (0) has been fixed. 

Here we consider, as usual, that the walker is initially located at the origin 
of the lattice, i.e., a„(0) = Vn except for n — 0. Then in the continuous limit 
we take 

A ± (C, 0) = oo(0)G (0 ± ± a.^G^), (53) 

where 



Gm(£) = AAexp 



4iu 2 



(54) 



with A/" a normalization constant that will be omitted in the following. As in 
[21] we are assuming that Eq. I|4(J|) is correct only for the long-wavelength com- 
ponents by taking an initial condition that "smears out" the lower-wavelength 
components. 

Now, by using Eqs. 1)46(1 and 1)53(1 . one easily obtains 



B ± {k,0) 



ao(0) ± o m (l) exp (mik) 



m=±l 

and with this, the result of (|4*5|) reads 



(55) 



B ± &r) 



aoW^ti^rli J] a m (l)Z ± [±(e-m),r] 



m=±l 



■>Ti\/p/<t>a 



where the functions Z ± (^',r) are 

r+oo 



(igexp 



iaq - -(3q 3 - (1 ± z 7 )<? 2 



7 



— sin(0 o r), 

VP 



2</> w 3 



sin(0 o r), 



cos{4>qt) - 1]. 



Their solutions read [231I2T] 



(56) 



(57) 



1 - a/3 - 7 2 ± 2i 7 



|/3|V3 

2 - 3a/3 - 6 7 2 3a/? + 2 7 2 - 6 



3/3 2 



3/3 2 



(58) 
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where Ai{z) is the Airy function. 

Finally, by using Eqs. I|56ll and (|46|l . we can write down the solution for the 
fields U^r) and D ± (£,t) 

L^K.t) = \u (0)Z ± {±Z,t}± ^(l)^ ± [±(C-m),r]}e ±2i iS9) 
I m=±l J 

D±(e,r) = JrfoW^Iie.rli 2 ^(l^IiCe-m)^]} 

I m=±l J 



The total probability of finding the walker in the position £ at time r can be 
easily calculated as 

V%t) = V u ^,r)+V d ^,r) (61) 
V u (t,T) = |«(e,r)| 2 = |u-+(t) + (-l)*l7-(i)| 2 , 

7> d (£,r) = |d(€,r)] 2 = |Z>+ (*) + (-1)*^ (t)| 2 - 
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